Draft version February 1 , 2008 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



DISCOVERY OF PULSATIONS AND A POSSIBLE SPECTRAL FEATURE IN THE 
X-RAY EMISSION FROM ROTATING RADIO TRANSIENT J1819-1458 

M. A. McLaughlin,'-' N. Rea,'*'' B. M. Gaensler," S. Chatterjee,'* F. Camilo,^ 
M. Kramer,^ D. R. Lorimer,'-^ A. G. Lyne,* G. L. Israel,' & A. Possenti* 

Draft version February 1, 2008 



o 
o 
o 

< 
OO 



o 



> 



OO 

o 

O 



ABSTRACT 

PSR J1819-1458 is a rotating radio transient (RRAT) source with an inferred surface dipole magnetic field 
strength of 5 x 10'^ G and a 4.26-s spin period. We present XMM-Newton observations of the X-ray counterpart 
of this source, CXOU J181939. 1-145804, in which we identify pulsations and a possible spectral feature. The 
X-ray pulsations are at the period predicted by the radio ephemeris, providing an unambiguous identification 
with the radio source and confirmation of its neutron star nature. The X-ray pulse has a 0.3-5 keV pulsed 
fraction of 34% and is aligned with the expected phase of the radio pulse. The X-ray spectrum is fit well by an 
absorbed blackbody with kT = 0. 14 keV with the addition of an absorption feature at 1 keV, with total absorbed 
flux of 1 .5 X 10"'^ ergs cm"^ s"' (0.3-5 keV). This absorption feature is well modeled by a Gaussian or resonant 
cyclotron scattering model, but its significance is dependent on the choice of continuum model. We find no 
evidence for any X-ray bursts or aperiodic variability on timescales of 6 ms to the duration of the observation 
and can place the most stringent limit to date of < 3 x 10"^ ergs cm"^ s"' on the absorbed 0.3-5 keV flux of 
any bursts. 

Subject headings: pulsars: individual (J1819-1458) — radio continuum: stars — stars: flare, neutron — 
X-rays: stars 



1. INTRODUCTION 

In 2006 February, a new class of neutron stars, the "Rotat- 
ing RAdio Transients" (RRATs) was reported (McLaughlin 
et al. 2006). These 1 1 objects, characterized by repeated dis- 
persed radio bursts, have periods ranging from 0.7 to 7 sec- 
onds and are located in the Galactic plane at 2 - 7 kpc dis- 
tances. Their periods are longer than those of most normal 
radio pulsars and similar to those of the populations of X- 
ray dim isolated neutron stars (XDINS; see Haberl 2007 for 
a review) and magnetars (see Woods & Thompson 2006 for 
a review). For the three sources with the highest pulse detec- 
tion rates, period derivatives, P, have been determined. If the 
measured P values are interpreted as due to magnetic dipole 
spin-down, they imply characteristic ages and magnetic field 
strengths in the general range of the normal pulsar population. 

There have been several suggestions put forward on the na- 
ture of this new class of neutron star One obvious suggestion 
is that the RRATs are related to pulsars which emit "giant 
pulses" (e.g. Knight et al. 2006). However, the RRATs with 
measured period derivatives do not appear to have high values 
of either magnetic field strength at the light cylinder or spin- 
down luminosity, both suggested as predictors of giant-pulse 
activity (Cognai'd et al. 1996; Knight et al. 2006). Zhang et 
al. (2006) suggest that the RRATs may be neutron stars near 
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the radio "death line" or may be related to "nulling" radio 
pulsars. However, the period derivatives measured for three 
RRATs do not place them near canonical pulsar "death lines" 
(e.g. Chen & Ruderman 1993) and, unlike most nulling pul- 
sars (e.g. Wang et al. 2007), we typically do not see more 
than one pulse from the RRATs in succession. Another in- 
triguing possibility is that the sporadicity of the RRATs is 
due to the presence of a circumstellar asteroid belt (Cordes 
& Shannon 2006; Li 2006) or a radiation belt such as seen 
in planetary magnetospheres (Luo & Melrose 2007). Or, per- 
haps, they are transient X-ray magnetars, a particularly rel- 
evant suggestion given the recent detection by Camilo et al. 
(2006) of transient radio pulsations from the anomalous X- 
ray pulsar XTE J1810-197. A final possibility is that they are 
similar objects to PSR B0656H-14, one of three middle-aged 
pulsars (i.e. "The Three Musketeers"; Becker & Truemper 
1997) from which pulsed high-energy emission has been de- 
tected (e.g. DeLuca et al. 2005). Weltevrede et al. (2006) 
convincingly show that if PSR B0656H-14 were more distant 
its emission properties would appear similar to those of the 
RRATs. Determining the reason for the unusual emission of 
the RRATs is important as population analyses show that their 
population may be up to several times greater than that of the 
normal radio pulsars (McLaughlin et al. 2006). Popov et al. 
(2006) show that the inferred birthrate of RRATs is consistent 
with that of XDINS but not with magnetars. 

While the radio emission properties of J1819-1458 are 
quite different from those of "normal" pulsars, it appears to 
be a rotating neutron star from which we detect radio pulses, 
and we henceforth give it the prefix "PSR". PSR J1819-1458 
shows the brightest and most frequent radio bursts of any of 
the RRAT sources. It has a 4.26-s period, relatively high in- 
ferred characteristic surface dipole magnetic field strength of 

5 X lO'-' G, characteristic age of 1 17 kyr, and spin-down lumi- 
nosity of 3 X 10^^ ergs s"' . The distance inferred from its dis- 
persion measure (DM) of 196±3 pc cm"^ is 3.6 kpc (Cordes 

6 Lazio 2002), with considerable (at least 25%) uncertainty. 
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PSR J1819-1458 is characterized by radio bursts at 1.4 GHz 
of average duration 3 ms, with bursts arriving randomly with 
a mean rate of one every ~ 3 minutes. X-ray emission was 
detected from this source in a serendipitous 30-ks Chandra 
ACIS-I observation toward the Galactic supernova remnant 
G15.9+0.2 (Reynolds et al. 2006)'. They found that the 
spectrum was well described by an absorbed blackbody with 
neutral hydrogen column density Nh = x 10^' cm~^ and 
temperature kT = 0.12±0.04keV, with an absorbed flux of 
~ 1 X 10"^^ ergscm"^ s"' between 0.3 and 5 keV. These prop- 
erties are consistent with emission from a cooling neutron star 
of age 10"*- 10^ years, broadly consistent with the charac- 
teristic age of PSR J18 19-1458. No evidence for bursts or 
variability was found. The time resolution of these Chandra 
observations was not sufficient to aUow a robust search for 
X-ray pulsations. 

We were awarded 43 ks of XMM-Newton time to further 
characterize the spectrum and search for pulsations. We re- 
port here on the results of these observations, in particular the 
detection of X-ray pulsations and a possible feature in the X- 
ray spectrum. This is the first detection of X-ray pulsations 
from any of the RRAT sources. In Section 2 we describe the 
observations and both the timing and spectral analyses. In 
Section 3 we discuss possible interpretations of our results, 
and present our conclusions in Section 4. 

2. OBSERVATIONS AND ANALYSIS 

The XMM-Newton observations were performed on 2006 
5 April, with a total observation time of 43 ks. The Euro- 
pean Photon Imaging Camera (EPIC) PN and MOS instru- 
ments were operated with medium filters and in Small Win- 
dow mode, providing a time resolution of 6 and 300 ms and 
effective livetimes of 71% and 97.5%, respectively. The data 
were reduced using the XMM-Newton Science Analysis Sys- 
tem (SAS version 7.0.0) and the most recent calibration files. 
Data from both the PN and MOS instruments were used for 
the timing and spectral analyses. We also also analyzed the 
Reflection Grating Spectrometer data but our target was too 
faint to be reliably detected. 

We filtered the observation for background flares, result- 
ing in effective on-source exposure times of 32.5 and 37.8 ks 
(23.0 and 36.9 ks including deadtime) for the PN and MOS in- 
struments, respectively. We detect a point source with J2000 
coordinates: right ascension a = 18''19™34^ and decUnation 
S = -14°58'03" (4" error in each coordinate). This is con- 
sistent with the radio-timing-derived position and with the 
more accurate position for the X-ray counterpart published 
in Reynolds et al. (2006). There is no evidence for extended 
emission, with the source brightness falling off as expected 
given the XMM-Newton point spread function. 

For both PN and MOS data, we extracted the source pho- 
tons within a 20" circular radius centered on the source posi- 
tion, which ensured extraction of more than 90% of the source 
counts. Following standard practice, the size of the extrac- 
tion region has not been corrected for in our final quoted 
fluxes and luminosities. The background counts were ex- 
tracted from four 20" circular regions centered on off-source 
positions (different for the PN and MOS instruments) in the 
same central CCD. 

2.1. Timing Analysis 

' As discussed by Reynolds et al. (2006), an association between the rem- 
nant and PSR J1819-1458 is extremely unlikely. 



For the timing analysis, we used all photons from PN and 
MOSl and MOS2 instruments with PATTERN < 12 (i.e. al- 
lowing for single, double, triple and quadruple events^"). This 
resulted in 2200 ±58 and 945 ±15 source and background 
PN counts and 1380 ± 41 and 230 ± 8 source and background 
MOS 1 plus MOS2 counts (normalized for a 20" circular ex- 
traction region). Although we can measure a period of 4.26 s 
from the arrival times of the radio bursts, the periodicity is not 
detectable in a Fast Fourier Transform of the radio time series 
due to the sporadic nature of the source. We therefore per- 
formed a periodicity search of all PN and MOS photons (i.e. 
with energies 0.2-15 keV) as a confirmation of the method 
through which the radio period was derived (see McLaughlin 
et ill. 2006). The arrival times were converted to the solar sys- 
tem barycenter and the test (Buccheri et al. 1983) was ap- 
plied. The most significant signal is detected with Zj = 144.5 
at a frequency of 234.564 ± 0.007 mHz (all errors are quoted 
at the 1(7 confidence level). Allowing for 4.3 x 10^ indepen- 
dent frequencies searched in the range 0.1-100 Hz, the prob- 
ability of chance occurrence of this signal in the absence of 
a real pulsation would be 1.8 x 10"^^, showing that the peri- 
odicity from this source would be easily detectable in a blind 
search of the X-ray photons. 

The ephemeris obtained through continued radio timing of 
PSR J1819-1458 wifli flie Parkes radio telescope at 1.4 GHz 
using the TEMPO software package predicts a barycentric 
frequency at the center of the observation (MJD 53830.87029) 
of 234.566244±0.000001 mHz, consistent with the measured 
X-ray frequency. Accounting for spin-down during the obser- 
vation is unimportant as the change in frequency due to the 
-3.16 x 10"^* Hz s"' frequency derivative over the 12-hr ob- 
servation is much less than the size of an individual frequency 
bin. The X-ray detection of periodicity shows that the radio- 
derived period is indeed the true period, and not a smaller 
common factor of the radio arrival times (see McLaughUn et 
al. 2006). 

In Figure 1, we present the 0.3-5 keV background- 
corrected X-ray pulse profile formed by folding barycentered 
photons from the PN and both MOS instruments using the 
radio ephemeris. The same good time interval file was used 
for the PN and both MOS instruments and backgrounds were 
subtracted separately. The different exposure times for the 
two instruments were accounted for. We chose to restrict this 
analysis to 0.3-5 keV as below 0.3 keV, the PN and MOS 
detectors are not weU calibrated and, above 5 keV, the back- 
ground dominates. In Figure 1, we also present the radio pro- 
file formed from adding individual bright radio pulses with 
the radio ephemeris. The X-ray pulse profile has a 0.3-5 keV 
background-corrected pulsed fraction of 34±6%, defined as 
iFmwc-Fn,in)/iF„ax+Fmm), where F„ca and F^in are flie mini- 
mum and maximum values of the X-ray pulse profile. It can 
be well modeled as a single sinusoid with xi ^ 1 .2 (17 d.o.f.). 
We find no significant dependence of pulsed fraction on en- 
ergy, with a pulsed fraction at energies 0.3- 1 keV of 28±7% 
and a pulsed fraction at energies from 1-5 keV of 49±10%. 

The peak of the X-ray profile, calculated by fitting a sinu- 
soid, is at phase 0.49 ± 0.04 and is aligned to the peak of the 
radio profile at phase 0.5. The 3 pc cm"^ uncertainty in the 
DM of 196 pc cm~^ amounts to a radio arrival time uncer- 
tainty at 1.4 GHz of only 13 ms, or 0.3% of the 4.26 s pulse 
period. 

Our PN data can be used to place limits on the existence 
http://xmm.viIspa.esa.es/extemal/xmm_user_support/documentation/uhb 
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of X-ray bursts from the source. Binning the data on various 
timescales to search for bursts of different widths we find no 
evidence for aperiodic variability on timescales ranging from 
6 ms to the duration of our observation. We used photons 
with energies up to 10 keV, important given the hard spec- 
tra of magnetar bursts (e.g. Gavriil et al 2004). Omitting 
frames with high particle background, we find that none of 
the 4x10^ frames contains more than two counts within the 
20" region, with three frames containing exactly two counts. 
These numbers are entirely consistent with Poisson statistics. 
The number of events in a single frame which would deviate 
from a steady flux at the 3-(t level is 3 photons. We there- 
fore adopt this as the upper limit on the X-ray flux of any 
burst of width 6 ms or less. Assuming a spectrum of the same 
shape as that fitted for the overall source (colunon 2 of Ta- 
ble 1 ; see Section 2.2) limits the observed fluence of any burst 
to < 8 X 10"'^ ergs cm"^ (0.3-5 keV), corresponding to an ab- 
sorbed flux Hmit of < 3 X 10"' ergs cm~^ s"' (0.3-5 keV) if 
we assume that the X-ray burst would last for 3 ms, the width 
of the radio pulse. This is a factor of five better than the limits 
placed by Reynolds et al. (2006). At a distance of 3.6 kpc, 
these limits correspond to a total energy of 2 x 10^"* ergs and 
luminosity of 1 x 10^^ ergs s~^ Therefore, any bursts, at least 
during the 12-hour duration of our observation, must contain 
much less energy than typical X-ray magnetar bursts (e.g. 
Woods & Thompson 2006) if their spectra is indeed similar 
to that described in colunm 2 of Table 1 . 

2.2. Spectral Analysis 

For the spectral analysis, we used PN photons with PAT- 
TERN < 4 (i.e. single and double events) and MOSl and 
MOS2 photons with PATTERN < 12. Source and background 
spectra were extracted from the same regions used for the tim- 
ing analysis and the spectral response matrices were created 
with the SAS mkrmf and mkarf tools, using the bad-pixel 
file built for our observation. The PN and MOS spectra were 
only used in the 0.4-2 keV energy range, a smaller range than 
that used for the timing analysis due to the greater dependence 
of spectral fitting on background spectra. In the spectral anal- 
ysis we used both rebinned and unbinned spectra. Spectra 
were rebinned for MOS and PN by a factor of two so that the 
energy resolution was not oversampled by more than a factor 
of three (using the specific response matrix built for this ob- 
servation). Furthermore, we rebinned in order to have at least 
30 counts per bin so that we could use the statistic" . 

We first modeled the PN spectrum. We tried several dif- 
ferent models and found that a single component fit was not 
possible. Fitting with a single absorbed'^ blackbody (as in 
Reynolds et al. 2006), we found Nh ~ 4.0 x 10^' cm'^ and 
kT ^ 0.14keV with a x^, > 2.0 (see Figure 2), while for an 
absorbed powerlaw, we obtained Nh ~ 1.3 x 10^^ cm"^ and 
r ~ 8.2 wifli a xl >21 (both fits have 53 d.o.f). These 
values are consistent with the Reynolds et al. (2006) best- 
fit parameters (0.5-8 keV) of Nh = TJ^ x 10^^ cm'^ and 
= 0. 12 ± 0.04 keV (for the absorbed blackbody model) and 
r ~ 9.5 (for the power-law model). An inspection of the resid- 
uals revealed that our high values were due to the pres- 
ence of strong spectral features around 1 keV and a weaker 
one around 0.5 keV. 

We carefully checked whether these features might be due 

' ' http://heasarc.nasa.gov/xanadu/xspec/manual/manual.html 

If not otherwise specified, abundances were assumed to be solar and 
fixed at the values in Lodders (2003). 



to calibration issues, to our source and background extraction 
regions or to residual particle flares and/or particles hitting the 
detector. We reliably excluded all of these issues by studying 
in detail the XMM-Newton calibration lines (see Figure 3), 
extracting the source and background photons from several 
different regions and investigating the spectrum using only 
PATTERN = counts (i.e. only isolated events). However, 
we note that the 0.5 keV feature is very close to the Oxygen 
edge energy, and an overabundance of Oxygen in the direc- 
tion of the source could be responsible. To further investigate 
this possibility, we fit our spectrum with an absorbed black- 
body model using three different photoelectric cross-sections, 
those of Balucihska-Church & McCammon (1992, 1998) and 
of Verner et. al. (1996). This was aimed at studying the 
dependence of the significance of our 0.5 keV line on the cho- 
sen photoelectric cross-section, which drives the shape of the 
edge. We found that the residuals around 0.5 keV remain in 
all cases, although less significantly if we use the Verner et. 
al. (1996) cross-section. For that reason, we assumed this 
cross-section for all our spectral fitting, and we onaitted all 
photons in the 0.50-0.53 keV energy range (two bins in our 
rebinned PN spectrum) from our modeling. We tentatively 
conclude that this line is due to the Oxygen edge, but only 
future deeper observations will unambiguously confirm this. 

Keeping the number of components as low as possible, 
we tried to model our rebinned PN spectrum with either an 
absorbed blackbody (XSPEC model bbody) or power-law 
(powerlaw), in addition to modeling the ~lkeV feature 
in several ways: leaving free the abundances of the most 
abundant elements in the interstellar medium (ISM) with 
lines around IkeV (e.g. Ne, N, Mg, Fe), adding a Gaus- 
sian function (XSPEC model gabs), a Lorentzian function 
(lorentz), an absorption edge (edge) or a cyclotron reso- 
nant scattering model (cyclabs). Among our trials, only the 
models reported in Table 1 gave satisfactory values of x?, ~ 1 • 
We show the blackbody plus Gaussian fit and residuals in Fig- 
ure 4. Note that there is evidence in the residuals for narrow 
lines within our broad 1 keV line, implying that the line may 
be a blending of narrower lines. However, the addition of 
several single narrow lines is not statistically significant, and 
future data with better statistics are necessary for investigating 
this issue. 

We ran Monte Carlo simulations of 2x lO'* spectra (see Rea 
et al. (2005, 2007) for details) to estimate the significance of 
the spectral feature (as suggested by Protassov et al. 2002)'^. 
The spectra were simulated using models corresponding to 
columns 2 and 4 from Table 1 with the absorbed blackbody 
parameters varying within their 3a errors (see Table 1), ab- 
sorption line parameters completely free to vary, and fixing 
the number of counts of each spectra at the PN count value 
for PSR J1819-1458. After having generated lO'' spectra for 
each of the two models, we counted how many of these sim- 
ulated spectra showed an absorption feature at any energy, 
width or depth, only due to statistical fluctuations. None of 
the simulated spectra presented a Gaussian absorption line at 
any energy or width with a depth tg > 30, and no cyclotron 
component was detected with a depth > 0.3 for any energy or 
width. We therefore infer the significance of the line at 1 keV 
to be > 1/10^, corresponding to a > 99.99% probability (i.e. 
at least 4a) that this line is not due to statistical fluctuations. 

Even though the F-test is still widely used to assess the significance of 
features in X-ray spectra, it has been shown that this is not statistically correct 
(Protassov et al. 2002). 
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TABLE 1 

Spectral fits for PSR J18 19-1458 with EPIC-PN 



Blackbody (BB) plus Neon BB plus Gaussian BB plus edge BB plus cyclotron 



Nh 

NNe 


59+0.06 
"•-'='-0.04 

6±1 


Eg 


75+0.12 
'-"-0.09 

1 1 i+fl.04 
^■^^-0.03 




.57±0.06 

"•^^-0.01 




"■'*i-0.08 

99+0.03 

"■^^-0.02 






tg 


21+0.03 
"■^^-0.06 

150±60 




0.67±0.14 


Wcy 


"•■"-0.06 

1.2±0.2 


kT 

Abs. Flux 
Unab. Hux 


(1 ,44+0.008 
"■^^^-0.006 

1 5"^-^ 

^•^-0.5 

15±2 




nfi"*'"^ 

U.lJU_QQQg 

1 5+0.6 
^•-'-0.8 

26^2 




150+''"''5 
1 5+°-3 




n ,44+0.008 
"■^^^-0.006 

1 5+0-3 

^■-'-0.4 

34!^ 


Xl (d.o.f.) 


1.20 (50) 




1.19(48) 




1.17(49) 




1.13 (48) 



Note. — Fluxes are calculated in the 0.3-5 keV energy range, and reported in units of 10"" ergs s"' cm"^. Nu is in units of lO^^cm"^ and N^e is in solar 
units (always assuming solar abundances from Lodders 2003). The photoelectric cross-section of Verner et al. (1998) has been used for all fits. The values of kT 
(blackbody temperature), Eg (Gaussian line energy), ctg (Gaussian line width), Ee (edge threshold energy), E^y (cyclotron line energy) and Wcy (cyclotron line 
width) are in units of keV. The Gaussian line depth tg, edge depth Te and fundamental cyclotron line depth dcy are dimensionless. Errors are at the la confidence 
level. XSPEC models used are (from left to right): vphabs*bbody, phabs*gabs*bbody, phabs*edge*bbody and phabs*cyclabs*bbody. 



To check whether the line was also present in the MOS data, 
we simultaneously fit the PN and both MOSl and MOS2 re- 
binned spectra (see Figure 5). All models in Table 1 fit the 
PN, MOSl and MOS2 data well. However, because of the 
lower number of MOS counts, our 1 ke V line was not signifi- 
cant in the MOS spectra, especially because the rebiiming we 
did in order to use the statistic left only a handful of bins 
within the line. As a further check, we then used the unbinned 
PN and MOS spectra and the C-statistic (Cash 1979) to as- 
sess the goodness of an absorbed blackbody fit, without the 
inclusion of any fine. We found that for Nh = 0.5 x 10^^ cm~^ 
and kT = 0.14keV, the C-statistic is 1009.776 using 969 PHA 
bins. A Monte Carlo simulation of the fit showed that the 
probabihty of having a C-statistic < 1010 is 98.0%. Note 
that only if the C-statistic probability is smaller then 50% can 
the model be accepted. This shows that an absorbed black- 
body model alone can not reproduce the data, and hence that 
a line is also present in the MOS spectra. We show in Figure 5 
the PN, MOS 1 and MOS2 spectra modeled with an absorbed 
blackbody plus a Gaussian (column 2 in Table 1), with a more 
severe PN rebinning (with at least 40 counts per bin), to show 
the agreement between the three instruments. Note that all the 
models in Table 1 satisfactorily fit the previous Chandra data 
(Reynolds et al. 2006), although the absorption line is not 
statistically significant in those data because of the reduced 
number of counts. 

Furthermore, we checked whether the spectral feature pa- 
rameters and significance were dependent on the spectral bin- 
ning or on the continuum model. We found that the spectral 
feature does not vary significantly when we vary the amount 
of spectral binning (see also Figure 5), but that it is rather de- 
pendent on the chosen continuum (as is generally the case 
for every spectral feature). In fact, if instead of a blackbody 
model, we assume a power-law continuum, the feature is still 
detectable with similar spectral properties but at a lower sig- 
nificance of 98.9% (roughly 2.5cr). 

Hence, if we assume a blackbody continuum model (rea- 
sonable for an isolated neutron star) we can show with high 
confidence that a spectral feature is present in our data. How- 
ever, our statistics do not allow us to statistically prefer one 
model to the others and reliably ascertain its nature. Further- 
more, if future, more sensitive observations reveal a different 
source continuum, the significance of the feature may need to 



be revised. 

We also performed phase-resolved spectroscopy over two 
intervals of 0.5 in phase centered on the peak and the mini- 
mum of the pulse profile, respectively. We did not find any 
significant spectral variability between these two phase inter- 
vals at the 2a level. However the low number of counts in 
each spectrum does not allow us to reach any reUable con- 
clusions about the pulse-phase dependence of the spectral pa- 
rameters. 

3. DISCUSSION 
3.1. Interpretation of the spectral feature 

If we accept that the spectral feature is significant, the two 
main possible interpretations are as an atomic line or as a cy- 
clotron resonant scattering line. 

If the fine is atomic, then it could be due either to the neu- 
tron star atmosphere or, but less probably, to a peculiar abun- 
dances in the ISM in that direction (or perhaps even a com- 
bination of the two). At 1 keV, it might then be due to the N 
edge (although improbable given its low abundance), Ne edge 
(around 0.9 keV), Ne XI K-edge (at 1.19keV), or some Fe-L 
Unes, including Fe XX, Fe XXI and Fe XXII or Fe XXIII (all 
around 0.9 keV). Note that Fe lines are common in other kinds 
of neutron stars, though mostly coming from accretion disks 
in binary systems rather than from the neutron stars them- 
selves (see e.g. Cottam et al. 2002). The structure present in 
the residuals (see Figure 4 middle panel), which in our data 
is not significant, might be due to a blending of narrow lines 
which we could only model with a broad line due to our lim- 
ited number of counts. 

If the feature is due to proton cyclotron resonant 
scattering, the magnetic field inferred would be Bey = 

1.6£',v(keV)/jc 10''* G, where yc = (1 -2GM /c^Ry/^ is the 
gravitational redshift factor (M and R are the neutron star 
mass and radius, respectively). Assuming canonical values 
for M and of 1 .4 Mq and 10 km, we find = 2 x lO'"* G. 
This is slightly higher than the dipolar surface magnetic field 
inferred through radio timing through the standard formula 
B=3.2x 10''^PP'/2 = 5 X 10^3 G. However, this expression 
can be considered an order-of-magnitude estimate as it as- 
sumes a purely dipolar field, a neutron star radius of 10 km, 
a moment of inertia of 10"*^ g cm^ and an angle between the 
spin and magnetic axes a = 90 degrees. Accounting for this 



DISCOVERY OF PULSATIONS AND A POSSIBLE SPECTRAL FEATURE IN THE X-RAY EMISSION FROM RRAT Jl 8 19-14585 



and assuming a = 30 degrees would make the timing- and 
cyclotron-inferred fields consistent. In addition, the width 
and depth of the Une are consistent with the predictions of 
Zane at al. (2001) for proton-cyclotron absorption in highly 
magnetized neutron stars. We cannot detect any harmonics to 
the fundamental cyclotron line because our spectrum is back- 
ground dominated above 2 keV. However, as observed in some 
accreting sources (Heindl et al. 2004) and investigated for iso- 
lated neutron stars by Liu et al. (2006), in some cases the first 
harmonic is as deep as the fundamental cyclotron frequency. 
Therefore, it is possible, although unlikely, that the 1 keV fea- 
ture is the first harmonic, with the 0.5 keV fundamental coin- 
cident with the depression in the spectrum that we have inter- 
preted as due to an overabundance of Oxygen. We note that 
the observed lines are unlikely to be due to electron-cyclotron 
absorption as the inferred magnetic field would be a factor of 
2000 lower (i.e. the ratio of the proton to electron mass) and 
incompatible with that measured through radio timing. 

3.2. Relationship to other classes of neutron stars 

Our detection of periodicity at the radio period shows that 
the X-ray source reported by Reynolds et al. (2006) is un- 
doubtedly the counterpart to PSR J1819-1458. The pulsed 
fraction and sinusoidal pulse shape are similar to what is 
observed for other middle-aged X-ray-detected radio pulsars 
such as B0656+14 (e.g. De Luca et al. 2005), which has been 
observed to have similar radio properties to the RRATs (Wel- 
tevrede et al. 2006). 

The thermal emission from PSR J1819-1458 is consistent 
with a cooling neutron star. However, the temperature from 
our blackbody fit appears slightly higher than temperatures 
derived from blackbody fits for other neutron stars of similar 
ages (see discussion by Reynolds et al. 2006). Note that it is 
possible that PSR Jl 8 19-1458 was born spinning at a sizable 
fraction of its present period of 4.26 s. In this case, as dis- 
cussed by Reynolds et al. (2006), the characteristic age could 
be a considerable overestimate and the inferred temperature 
could be completely consistent with its age. Note that char- 
acteristic ages have been shown to be misleading for several 
other pulsars (e.g. Gaensler & Frail 2000; Kramer et al. 2003). 

Including PSR J1819-1458, eight high-magnetic field ra- 
dio pulsars (i.e. B > 1 x lO'"' G) have now been observed at 
X-ray energies. Two, J1846-0258 and B1509-58, are bright 
non-thermal sources (Mereghetti et al. 2002; Cusumano et al. 
2001), as expected given their young ages (less than 2000 yr). 
PSR 111 19-6127 is a bright thermal X-ray emitter with un- 
usual properties including a large pulsed fraction and narrow 
pulse (Gonzalez et al. 2005). PSR J1718-3718, with mag- 
netic field of 7 X lO'^ G, has been detected at X-ray energies, 
but the faintness of the counterpart does not allow detailed 
spectral modeling or a constraining limit on pulsed fraction 
(Kaspi & McLaughlin 2005). No X-ray emission has been 
detected from PSRs J1814-1744, B0154H-61 or J1847-0130, 
which has the highest inferred surface dipole magnetic field 
(9 X 10'^ G) measured to date for any radio pulsar (Pivovaroff 
et al. 2000; Gonzalez et al. 2004; McLaughlin et al. 2003). 
Radio pulsar X-ray emission properties clearly vary widely, 
even for objects with very similar spin-down properties. 

While the spectrum and luminosity of PSR J 1 8 19-1458 ar- 
gue against a relationship with magnetars, it is possible that 
PSR J1819-1458 could be a transition object between the pul- 
sar and magnetar source classes. The soft X-ray spectrum 
does have a comparable temperature to the quiescent state of 
XTE J1810-197 {kT ~ 0.15-0.18 keV; Ibrahim et al. 2004; 



Gotthelf et al. 2004). However, the radio emission character- 
istics of these two neutron stars are quite different. 

While resonant cyclotron features are regularly observed 
from X-ray binary systems (e.g. Truemper et al. 1978; Naka- 
jima et al. 2006), the detection of such features from iso- 
lated neutron stars is quite unusual. Sanwal et al. (2002) & 
Mereghetti et al. (2002) discovered harmonically spaced ab- 
sorption lines from IE 1207.4-5209, a radio-quiet X-ray pul- 
sar with a 424-ms spin period and timing-derived characteris- 
tic age and inferred surface dipole magnetic field strength of 
3 X 10^ yr and 3 x 10^^ G, respectively. Analyzing a deeper 
observation of the source, Bignami et al. (2003) attributed 
these fines to electron-cyclotron absorption. However, the sig- 
nificance of two of the lines present in its spectrum has been 
strongly questioned (Mori et al. 2005). 

Broad absorption lines, similar to those seen for 
PSR J1819-1458, have been observed for six out of seven 
XDINS (see reviews by van Kerkwijk & Kaplan 2007 and 
Haberl 2007). For most of these neutron stars, the lines can 
be interpreted as due to neutral hydrogen transitions in highly 
magnetized atmospheres. Ho et al. (2003) and Van Kerkwijk 
& Kaplan (2007) argue that the transition energy is similar to 
the proton cyclotron energy for magnetic fields of the order 
of PSR J1819-1458's. The X-ray spectrum of PSR J1819- 
1458 is very similar (although with a slightly hotter blackbody 
temperature) to the XDINSs, although so far no convincing 
evidence for radio bursts has been detected for any of those 
thermally emitting neutron stars (Kondratiev et al., in prepa- 
ration). 

One outstanding question is why absorption Unes of this 
kind, whether due to the atmosphere or to cyclotron resonant 
scattering, have been observed from only a handful of X-ray 
emitting isolated neutron stars. As suggested by Mereghetti 
et al. (2002), the age of the neutron star could be one key fac- 
tor. Young objects are dominated by non-thermal emission, 
but older ones may be too faint for X-rays to be detectable, 
making X-ray bright, middle-aged pulsars the best candidates 
(as is the case of the XDINSs and of PSR 1819-1458). Note, 
however, that no such absorption lines have been found for 
the X-ray bright, middle-aged PSRB0656+14, despite deep 
searches both with Chandra (Marshall & Schultz 2002) and 
XMM-Newton (De Luca et al. 2005) and the theoretical pre- 
dictions for cyclotron and/or atmospheric features expected in 
its emission. The explanation could well depend on the view- 
ing geometry. 

More sensitive observations are necessary to confirm 
the presence of absorption features in the spectrum of 
PSR J1819-1458. Longer X-ray observations are needed to 
understand whether the broad 1 keV line could be a blending 
of narrower features. This would then argue for an atomic 
origin for the line. Longer observations will also allow phase- 
resolved spectroscopy. Simulations show that a 100-ks obser- 
vation would allow us to achieve both of these goals. Because 
the strength of the cyclotron line depends on the angle be- 
tween the observer and the magnetic field, we expect phase 
variation of cyclotron features. If the feature we detect is in- 
deed due to proton-cyclotron absorption, it provides an invalu- 
able means of testing the assumptions implicit in the charac- 
teristic magnetic fields derived through radio timing measure- 
ments and a unique independent measurement of the magnetic 
field of an isolated neutron star. 

4. CONCLUSIONS 
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We have discovered X-ray pulsations at the 4.26-s period 
inferred from radio timing of PSR J1819-1458. The prop- 
erties of these pulsations are similar to those observed for 
other middle-aged radio pulsars detected at X-ray energies. 
While the RRATs are characterized by sporadic radio emis- 
sion, we do not detect any X-ray bursts or aperiodic varia- 
tions throughout the observation and can place the most strin- 
gent limit to date of < 3 x 10"^ ergs cm~^ s~' on the absorbed 
0.3-5 keV flux of any X-ray bursts. We have characterized 
the spectrum of this source and find that it is well-described 
by an absorbed blackbody with kT = 0.14 keV in addition to 
an absorption line around 1 keV, with total absorbed flux of 
1.5 X lO"'-' ergs cm"^ s~' (0.3-5 keV). We note, however, that 
the presence of this absorption feature is highly dependent on 
the choice of continuum model and needs further X-ray ob- 



servations to be confirmed. This object is the only RRAT so 
far to be detected at X-ray energies. X-ray observations of the 
other objects in this source class are essential for a complete 
picture of how they relate to other neutron star populations. 
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DISCOVERY OF PULSATIONS AND A POSSIBLE SPECTRAL FEATURE IN THE X-RAY EMISSION FROM RRAT Jl 8 19-14587 
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Fig. 1.— Top: X-ray pulse profile (0.3-5 keV) for 23-ks of EPIC-PN data and 37-ks of MOSl and MOS2 data on PSR J1819-1458. The dotted line shows 
the best-fit sinusoid and the vertical dashed line indicates the phase of the peak of the radio pulse. Bottom: radio profile formed from 114 pulses detected in 6 hr 
of observation at 1.4 GHz with the Parkes telescope in Australia. The non-uniform baseline is likely due to instrumental digitization effects. Both profiles have 
been folded using the radio ephemeris. In both plots, two cycles of pulse phase are shown. 
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Fig. 2. — Top: XMM-Newton EPlC-PN spectrum of PSR J1819-1458, modeled as an absorbed blackbody. The points indicate the data, while the solid line 
shows the corresponding best-fit model. Bottom: residuals of the model. 
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Fig. 3. — The effective area, or total response efficiency, for EPIC-PN with the medium filter versus incident photon energy. Note that there is no instrumental 
feature in the effective area curve at 1 keV. The calibration feature near 0.5 keV is too weak to explain the feature in our spectrum close to this energy (see Haberl 
et al. 2006 for more details on XMM-Newton calibration.) 



DISCOVERY OF PULSATIONS AND A POSSIBLE SPECTRAL FEATURE IN THE X-RAY EMISSION FROM RRAT Jl 8 19-14589 
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Fig. 4.— Top: XMM-Newton EPIC-PN uF^, spectram of PSR J1819-1458, modeled as an absorbed blackbody and Gaussian line (column 2 of Table 1). The 
points indicate the data, while the solid line shows the corresponding best fit model. Top line shows the absorbed blackbody model component of the fit alone. 
We have omitted all photons with energies from 0.5-0.53 keV as we believe that the feature in this range is due to the Oxygen edge (see Section 2.2). Middle: 
residuals of the absorbed blackbody and Gaussian line model. Bottom: residuals of the model without the inclusion of the Gaussian Une. 
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Fig. 5.— Top: XMM-Newton EPIC-PN (filled circles), MOSl (black cross), and M0S2 (triangles) spectra for PSR J1819-1458, modeled as an absorbed 
blackbody and Gaussian line (column 2 of Table 1). Grey lines show the fitted model. The PN data have been binned with 40 counts per bin and the MOSl and 
MOS2 data with 30 counts per bin to show that the best-fit model is not dependent on the spectral binning parameters. We have omitted all photons with energies 
from 0.5-0.53 keV as we believe that the feature in this range is due to the Oxygen edge (see Section 2.2). Bottom: residuals of the above fits. 



